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Abstract 



A molecular dynamics simulation is performed for a supercooled liquid of 
rigid diatomic molecules. The time-dependent self and collective density cor- 
relators of the molecular centers of mass are determined and compared with 
the predictions of the ideal mode coupling theory (MCT) for simple liquids. 
This is done in real as well as in momentum space. One of the main results 
is the existence of a unique transition temperature T c , where the dynamics 
crosses over from an ergodic to a quasi-nonergodic behavior. The value for T c 
agrees with that found earlier for the orientational dynamics within the error 
bars. In the /3-regime of MCT the factorization of space- and time depen- 
dence is satisfactorily fulfilled for both types of correlations. The first scaling 
law of ideal MCT holds in the von Schweidler regime, only since the validity 
of the critical law can not be confirmed, due to a strong interference with 
the microscopic dynamics. In this first scaling regime a consistent description 
within ideal MCT emerges only, if the next order correction to the asymptotic 
law is taken into account. This correction is almost negligible for q = q max , 
the position of the main peak in the static structure factor S(q), but becomes 
important for q = q m i n , the position of its first minimum. The second scal- 
ing law, i.e. the time-temperature superposition principle, holds reasonably 
well for the self and collective density correlators and different values for q. 
The a-relaxation times Tq and r q follow a power law in (T — T c ) over 2-3 
decades. The corresponding exponent 7 is weakly q-dependent and is around 
2.55. This value is in agreement with the one predicted by MCT from the 
value of the von Schweidler exponent but at variance with the corresponding 
exponent 7 ~ 1.6 obtained for the orientational correlators c[ s \t) and C\(t), 
studied earlier. 
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I. INTRODUCTION 



Although several interpretations of the glass transition exist (see e.g. Ref. [1]]), the only 
microscopic approach leading to a variety of predictions is based upon the mode coupling 
theory (MCT). This theory, which was proposed first by Bengtzelius, Gotze and Sjolander 
Oj and Leutheusser Q , provides an equation of motion for the normalized density correlator 



of a simple liquid. Note that F depends on q = \q\ only, due to the isotropy of the system. In 
its idealized version, MCT predicts the existence of a dynamical transition at a critical tem- 
perature T c (or a critical density n c ) from an ergodic to a nonergodic phase, corresponding 
to a liquid and a glass, respectively. The nonergodicity parameter (NEP) 

/(?)=lini%t) (2) 

t — >oo 

serves as an order parameter for that transition. f(q) may change either continuously (type- 
A transition) or discontinuously (type-B transition) at T c . For the structural glass transition 
only the latter is relevant. 

The test of the predictions of the ideal MCT has challenged both, experimental work and 
computer simulations. Most of the experiments have determined the intermediate scattering 
function F(q, t) or its time-Fourier transform, the dynamical structure factor (coherent and 
incoherent part). This yields information on the dynamics of the translational degrees of 
freedom (TDOF), only. A great part of these investigations was reviewed by Gotze and 
Sjogren || and in Refs. ||. 

Since most of the glass formers are molecular systems, it is also important to study 
the role of the orientational degrees of freedom (ODOF) and their coupling to the TDOF. 
A convenient way for this is, e.g., the use of dielectric spectroscopy. Recent results from 
dielectric measurements || are partly consistent with the predictions of MCT for simple 
liquids. The interpretation of light scattering experiments is not so obvious. Since TDOF 
and ODOF may both contribute to light scattering, it is not so easy to separate the orien- 
tational part, (see e.g. Ref. 0). Concerning computer simulations, only a few molecular 



liquids have been studied so far f8Hl5|. In a recent molecular dynamics [MD]-simulation [16 
the present authors investigated a liquid of diatomic, rigid molecules with Lennard- Jones 
interactions. Apart from linear molecules with head-tail symmetry, this is the simplest 
choice for a molecular system. Because dielectric and light scattering measurements yield 
essentially only information for q rs 0, we restricted ourselves to orientational correlators 
with q = and to the translational and rotational diffusion constant D and D r , respectively 
16fl . The main result we have found is the existence of a power law dependence of all the 



orientational relaxation times r/ (Z = 1 — 6), t\ (I — 1) and of D on temperature with a 
single critical temperature T c = 0.475, in agreement with ideal MCT. Slightly above T c the 
numerical results for r/ 5 '', tj and D deviate from a power law, due to ergodicity restoring 
processes. These processes, which are not considered by the ideal MCT, can be accounted 



for by its extended version [17|. The observation that the temperature dependence of D r 



does not fit into this scheme, but bifurcates from D at a temperature significantly above T c 
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and then obeys an Arrhenius law, rather than a power law, has been a further important 
result. 

The main purpose of the present and the following paper [I8| is the extension of our 
recent results [[U| to finite q and to perform a thorough test of the predictions of the ideal 
MCT. The present paper is restricted to the TDOF only, whereas the following paper [18 



considers orientational correlators, which explicitly contain the coupling between ODOF and 
TDOF. Although our simulation investigates a molecular system, we make a comparison 
with the predictions of MCT for simple liquids. The ideal MCT for simple liquids has 
been recently extended to a dumbbell molecule in a simple liquid |TjJ] and to molecular 
liquids of linear pO|j2T[] and arbitrary shaped molecules p2|. So far, the investigations of 



these molecular MCT-equations were restricted to the calculation of the NEP ||19| , |21| , |22 



For linear molecules without head-tail symmetry, it follows from these equations |21[ that 
the TDOF and ODOF freeze at a single critical temperature T c , 




consistent with our recent 



results for q 



m 



The investigation of the time dependent molecular MCT-equations, 
which would allow comparison with our MD-results, will be done in future. 

The outline of this paper is as follows. The next section will review the predictions of the 
ideal MCT (for simple liquids) which will be tested. In section |TTl| the model as well as some 
details of the computer simulation are discussed. In section [TV] we present our MD-results 
and the final section contains a discussion of these results and our main conclusions. 



II. MODE-COUPLING THEORY 

In this section we will give a short summary of those predictions of the ideal MCT 
(for simple liquids), which will be compared with our MD-results. For details the reader 
may consult the review papers For completeness we also would like to mention the 

MCT-approach by use of nonlinear fluctuating hydrodynamics p3| , p4 . 



Ideal MCT predicts the existence of a dynamical transition at T c from an ergodic to a 
nonergodic phase which is an ideal glass transition. For temperatures close to T c , MCT 
predicts the existence of two scaling laws for F(q, t) with time scales t a and r(^> t a ) , where 
r is the a-relaxation time. In the first scaling regime, i.e. for t <C t <C r, the density 
correlator takes the form: 

F(q,t) = r(q)+h(q)G(t) (3) 

with f c (q) the NEP at T c and h(q) the critical amplitude, to is a microscopic time of the 
order of 10~ 13 sec. The so-called /5-correlator (or critical correlator) G(t) obeys the first 
scaling law: 

G(t) = c a g ± (t/Q (4) 

with the correlation scale c CT = | cr | 1//2 and the separation parameter a = oq(T — T c ), where 

o"o > 0. The ex-independent master functions g±(& < 0) are solutions of a certain scaling 
equation. The corresponding /3-relaxation time scale t a is given by: 

Uoc \T-T C \~^ , T<T C (5) 
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For G(t) one obtains the following asymptotic power laws: 



The critical law, upper equation of (6), holds above and below T c , whereas the von Schweidler 
law, lower equation of (6), is valid only above the transition point. Below T c , G(t) decays 
to a constant for t ^> t a . Both exponents a and b are related to the exponent parameter 
A(0 < A < 1) by: 

r(i-q) 2 = = ni + bf 

T(l-2a) r(l + 26) ' {) 

with T the Gamma function. From (7) one gets that < a < ~ and < b < 1. A is 
determined by the static correlators at T c . The result (3) which states that in the first 
scaling regime the q- and t-dependence factorizes is one of the most important predictions 
of the ideal MCT. Due to this factorization, Eq. (3) can easily be transformed to real space: 

<f ) (r,t) = f c (r) + H(r)G(t) (8) 

where <p(r,t), f c (r) and H(r) are the Fourier transform of F(q,t), f°(q) and h(q), respec- 
tively. 

It is important to realize that the result (6) only holds asymptotically for T — ► T c . The 
next order correction to both of these asymptotic laws was recently derived and calculated 
for a system of hard spheres by Franosch et al. [ 25| . Corrections to the von Schweidler law 

m 

F(Q,t) = /%) - Kq)(t/r) b + h^{q){t/rf b - . . . (9) 



have already been studied for hard spheres by Fuchs et al. ||27|| . There it was demonstrated 
that these corrections may be important. The expansion (9) is valid for t a t r, where 
the a-relaxation time scale r is given by 

r(T) oc (T — Tc)" 1 , T>T C (10) 

with 

7 = — + — 11 

1 2a 2b y 1 

The translational diffusion constant D is predicted to scale like r _1 . Therefore it is: 

D{T) oc (T — T c y , T > T c . (12) 



Corrections to the critical law lead to [25] : 

F(q,t) = f c (q) + h(q)(t/t a )- a + U 2 \q){t/t a )- 2a + ... (13) 

In the second scaling regime, i.e. for t of the order of r, a master function F q (t) exists such 
that 
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F(q,t,T) = F q (t/r(T)) 



(14) 



In glass science the result (14) is called time-temperature superposition principle. The 
expansion of the r.h.s. of (14) with respect to t/r yields (9). Equation (14) represents the 
second scaling law. 

Finally we mention that all these MCT-results hold for the self part of the density 
correlator as well. 



III. MODEL AND DETAILS OF THE SIMULATION 

The model we investigate is a one-component system of rigid diatomic molecules. Each 
molecule is composed of two different Lennard- Jones particles, in the following denoted by 
A and B, which are separated by a distance d = 0.5 and each of which has the same mass m. 
The interaction between two molecules is given by the sum of the interaction between the four 
particles which is given by the Lennard- Jones potential V a p(r) = 4e Q , j g{(cx Q , / 3/r) 12 — (cr a/ 3/ r ) 6 } 
where a, /3 G {A,B}. The Lennard- Jones parameters are given by: oaa = <?ab = 1-0, 
a BB = 0.95, €aa = £ab = 1-0 and e B B = 0.8. In the following we will use reduced units and 
use oaa as the unit of length, €aa as the unit of energy (setting k B = 1) and (a\ A m/A'&eAA) 1 
as the unit of time. If the atoms are argon-like this time unit corresponds to approximately 
0.3 ps. 

In order to make the simulation more realistic we did it at constant external pressure 
p ext =1.0. The length of the equilibration runs always exceeded the typical relaxation time 
of the system at the temperature considered, which allows us to conclude that in the subse- 
quent production runs we were investigating the equilibrium dynamics of the system. The 
temperatures we investigated are T = 5.0, 3.0, 2.0, 1.4, 1.1, 0.85, 0.70, 0.632, 0.588, 0.549, 
0.520, 0.500, 0.489 and 0.477. The total number of molecules was 500 and, in order to 
improve the statistics of the results, we averaged at each temperature over at least eight 
independent runs. For more details see jL6 . 



IV. RESULTS 

For a clearer presentation of our results, this section is divided into two subsections, 
where the first is devoted to the static and the second to the dynamical properties. The 
latter is divided into two parts again, which are real space and g-space representation. 



A. Static properties 

First of all we remark that thermodynamic quantities like, e.g. the average density and 
the enthalpy do not exhibit any signature of a singular behavior for the full temperature 
range 0.477 < T < 5.0, we have investigated. From this we conclude that the observed 
slowing down is probably not related to an approach of the system to a critical point of a 
second order phase transition. 

The structural properties are one of the most interesting static features of a supercooled 
liquid. Figure 1 shows the static structure factor S(q) of the center of mass positions for 
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different temperatures as a function of q. Its g-dependence has the typical behavior expected 
for a liquid, with a main peak at q ma x — 6.5 and a first minimum at q m in — 8-15 for the lowest 
temperature. With increasing temperature the peak positions shift to smaller q values, due 
to an increase of the average distance between the molecular centers, and the peaks become 
broader. 

A comparison of S(q) for T = 0.477 with the partial structure factors S a p{q), where a, (3 
refers to atom A and atom B, is given in Fig. 2. Whereas the main peak in S(q) and S a p(q) 
is essentially at the same position, the remaining g-dependence is quite different for the 
various correlators. For the pre-peak of Saa(q) and Sbb(q) at Q ~ 3.1, which could indicate 
a "medium" range order, we have not found a conclusive interpretation. Concerning the 
microscopic structural properties it is more convenient to consider the atomic radial pair 
distribution functions g a p(^) and the corresponding function <7 COT (r) for the center of mass 
positions, which are shown in Figs. 3 and 4, respectively. There is no significant difference 
between the g a p{r) probably due to the fact that the choice of our model parameters (T^ 
does not introduce a strong asymmetry between particles of type A and particles of type B. 
Besides a first and a splitted second nearest neighbor peak at about r\ ~ 1.05 and r-i « 1.9, 
r' 2 ~ 2.1, respectively (Fig. 3), an intermediate peak occurs at r[ ~ 1.5. The peak at r' 2 is 
also found for g cm , but its main peak is at r'{ ~ 1.25 (cf. Fig. 4). The peaks in g a/3 at r[ and 
r 2 are missing for g cm . The splitted first nearest neighbor peak of g cm can be attributed to 
cross- and T-configurations, as also found for other molecular liquids (see e.g. |2^|). We do 
not dwell more on these results, since our main concern are the dynamical features in the 
strongly supercooled regime. The main reason we have presented S(q) is to compare its q- 
dependence with the static correlation functions of the orientational correlators presented in 
part II. We already anticipate that this g-dependence can be rather different for the various 
correlators. 



B. Dynamical properties: real space 

To start we discuss the self motion of a molecule. One of the important transport 
coefficient is the translational diffusion constant D which is obtained from the mean squared 
displacement by: 

1 N 

^fc^E((a)-4(0)f) , (15) 

where N is the number of molecules. The temperature-dependence of D is presented in Fig. 5 
on a double-logarithmic scale. The diffusion constant follows a power law, as predicted by 
ideal MCT [cf. Eq. (12)], over a surprisingly large range of four decades in D. As critical 
temperature we obtain T c = 0.475 and for the corresponding exponent jjj = 2.20. With use 
of Eqs. (7) and (11), this yields the von Schweidler exponent be = 0.69, the critical exponent 
cld = 0.34 and the exponent parameter = 0.67. At the two lowest temperatures the 
deviations of the numerical values for D from the power law, which has been mentioned in 
the Introduction, can clearly be seen. The inset of Fig. 5 depicts the time dependence of 
the mean squared displacement for all temperatures, from which D(T) was deduced. 

More information on the self motion is contained in the self part of the van Hove corre- 
lator: 
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1 N 



Gs(r, t) = (-Y, 5 ( r - - 4(0)1) ) (16) 



n=l 



The r-dependence of G s (r,t) is shown in Fig. 6 for times which are equidistant on a loga- 
rithmic time axis and for the lowest temperature. With time, the 5-peak at t = broadens 
and for t — > oo and r — > oo, i.e. in the hydrodynamic limit, it approaches a Gaussian 
distribution: 

In the time span 5.0 < t < 10 3 , where the mean squared displacement for the lowest 
temperature exhibits a plateau (cf. the inset of Fig. 5), the r-dependence of G s (r,t) varies 
only weakly with time, which indicates the validity of the factorization [Eq. (8)]. To check 
this, we have plotted in Fig. 7 (G a (r,t) — G s (r,t'))/(G s (r',t) — G s (r',t')) as a function of r 
for fixed r' and t' and various t-values. If Eq. (8) is valid, this ratio must be independent of 
t and t'. For the times t and t' we have chosen in Fig. 7, this is indeed reasonably fulfilled. 
That this is not true for the full time span, can be inferred from the inset of Fig. 7. 

Another interesting conclusion can be drawn from the r-dependence of G s for the largest 
time. As can be seen from Fig. 6, at no time there is any indication for the presence of a 
secondary peak at r ~ 1, the nearest neighbor distance. From the absence of such a peak it 
is usually concluded, see e.g. Ref. [p9fl , that no hopping processes are present. But for the 
present system we have strong evidence, see Fig. 5, that at the lowest temperature hopping 
processes are indeed present. Thus we conclude that hopping processes do not necessarily 
lead to a secondary peak in G s . 

The r-dependence of G s is not Gaussian for a large time regime. The deviation from a 



Gaussian can be quantified by the non-Gaussian parameters oc n (t), n = 2, 3, . . . [30]. Fig. 8 
shows 

a 2 (t) = ^%-l . (18) 

In case of a Gaussian process a 2 (i) vanishes. For t — > and t —>■ oo a 2 (t) goes to zero. 
But in between there exists a time regime for all temperatures where a 2 is substantially 
different from zero. The increasing part of a 2 seems to fall onto a master curve, as already 
found earlier fll3l^9"l . The time t max (T), where a 2 (t) obtains its maximum, obeys a power 
law (T — T c )~ la (see inset of Fig. 8) with 7 a pa 2.1 which coincides with jjj, derived from 



D(T), in agreement with the results of Sciortino et al. This suggests that a 2 (t) for 

t ~ t max is determined by the ct-relaxation as it follows from ideal MCT. Exactly the same 
conclusions also hold for a 3 (t). 

The collective dynamics can be studied by use of the distinct part of the van Hove 
correlator 

= (a73T £ 6(r - \x n (t) - S m (0)\)\ (19) 
which is shown in Fig. 9 for various times and again for the lowest temperature. Its r- 



and t-dependence is quite similar to that found for a binary liquid |[29|| . The test of the 
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factorization [Eq. (8)] can be done as for the self part. The corresponding ratio is presented 
in Fig. 10 for the same time regime as in Fig. 7. Again a reasonable collapse onto one master 
curve occurs, in agreement with ideal MCT. 



C. Dynamical properties: g-space 

Although the visualization of local structural properties favors the use of a real space 
representation of the correlators, most experimental results are obtained in g-space. In 
addition the translational invariance of the interactions also suggests a theoretical description 
by use of g-dependent correlators. For instance the memory part of the MCT-equation for 
the density correlator F(q, t) is "diagonal" in g-space but involves a convolution in real 
space. 

Let us start with the self part F s (q,t) of the collective density correlator F(q,t) [cf. 
Eq. (1)], which is presented in Fig. 11 for q = q max = 6.5 , the position of the main peak 
in S(q) (cf. Fig. 1). We remind the reader that q is measured in units of oaa-, the position 
of the Lennard- Jones minimum for the A-particles. The small bump in F s (q,t) at t s = 12, 
which is also present in F(q, t) (cf. Fig. 14 and 15), is due to a sound wave and t s is the time 
for the sound wave to traverse the box size of our sample. This effect was also observed by 



Lewis et al. |TIJ and has recently been shown to be much more pronounced in strong glass 
formers f31~ [. 

For the highest temperature, T = 5.0, the decay of F s (q,t) is essentially exponential. 
With decreasing temperature the relaxation crosses over into a two-step process, as predicted 
by MCT. For the lowest temperature, T = 0.477, a "quasi-plateau" can be seen. Its height 
is a measure of the critical NEP f^(q). The reader should note that F s (q max ,t) decays to 
zero for large times, indicating that the length of the runs is large enough to equilibrate the 
system. In order to test MCT, we have determined the a-relaxation times t^ s \T) from the 
condition F s (q, t' ! ') = 1/e. Its temperature dependence is shown in Fig. 12 (squares) for 
q = q m ax = 6.5 and q = q m i n = 8.15. Using T c = 0.475, as determined from D(T), both 
relaxation times follow a power law over about three decades. The corresponding 7- values are 
practically the same 7^ = 2.56, but differ significantly from 7^ = 2.20. Such a discrepancy 
between 7^ and 7^ has been observed already before and indicates that the MCT 



prediction, that the two exponents should be equal, is not valid for these systems. 7 ( - s ^ = 2.56 
yields the von Schweidler exponent b = 0.55 and the exponent parameter A = 0.76. We also 
note that a power law fit with T c as a free parameter yields a value for T c which coincides 
with T c = 0.475 to within 2%. 

The relaxation time t^(T) can now be used to rescale time in order to test the time- 
temperature superposition principle [Eq. (14)], which is done in Fig. lib. We find that 
the relaxation curves fall indeed onto one master curve for the lowest temperatures. The 
late stage relaxation can be fitted well by a Kohlrausch- Williams- Watts-law (KWW), i.e. 
F s (q,t) = Aexp -(t/r^f , dashed curve. For higher temperatures (cf. e.g. the curve for 
T = 5.0) no such scaling exists. We also note that such a scaling was not possible for certain 
types of orientational correlation functions [IB], which shows that this prediction of MCT 



is not a trivial one. 

Having demonstrated the validity of the second scaling law [Eq. (14)], we can test whether 
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the von Schweidler law, including the next order correction [Eq. (9)] fits the master curve 
well in the late /3-relaxation regime. As can be seen from Fig. lib this type of fit works 
very well (dotted line). In practice, this has been done for F s (q, t) at the lowest temperature 
by using b = 0.55 (deduced from 7^). The wave vector dependence of fg(q), h s (q) and 
h^ s (q) is shown in Fig. 13. Here a comment is in order. The fit of the data with the 
von Schweidler law including the next order correction yields r~ b h s (q) and T~ 2b h^ s (q) [cf. 
Eq. (9)]. Since the g-independent a-relaxation time scale r{T) can only be determined up to 
a T-independent factor, the same is true for h s (q) and h^ s (q). In Fig. 13 we therefore show 
h s (q) = r~ b h s (q) and h^ s (q) = r~ 2b h g 2 ^ s (q). The q- variation of these quantities is quite 
similar to that for hard spheres ||27|| . Particularly we also find a zero of h^ s (q) at a finite 
q-value go ~ 3.5 with h^ s (q) < 0(> 0) for q < q (q > q ). Taking into account that q for 
the hard sphere system is given in units of the diameter d HS = 1, we can deduce q$ ~ 12 
from Ref. ||27|| . That this is about three times larger than our value, is partly due to a larger 
effective diameter d e ff « 1.5 of our molecules, compared to d HS = 1. 

We now turn to the collective dynamics as obtained from the density correlator F(q,t), 
which is presented, respectively, for q = q max and q = q m i n in Figs. 14a and 14b. The 
time-dependence looks similar to that of F s (q,t), but the height of the "quasi-plateau" is 
much lower and the stretching of the relaxation is more pronounced for the case of q = q m in- 
The a-relaxation time r q (T) has been determined from the condition F(q,r q ) = 0.1, since 
using the e -1 definition of r would lead to an underestimation of the a-relaxation time 
for q = q m in, since the plateau is so low. (We note that choosing e -1 instead of 0.1, for 
the case where both definitions can be used, or even using a KWW-fit to deduce r q (T), 
leads to essentially the same temperature-dependence). The temperature dependence of 
r q is shown in Fig. 12. Taking again T c = 0.475 as given, T q (T) obeys a power law, with 
lq ma x = 2-57 ~ 7* s - ) and 7<j m4n = 2.47. It is remarkable that the range for the power law for 
q = qmin is almost one decade less than for q = q max , which shows that different correlators 
reach the asymptotic regime at different temperatures. 

Fig. 15 shows that also for these correlation functions the second scaling law holds and 
that it holds better for q max than for q m i n . The late stage relaxation can again be fitted 
well by a KWW-law (dashed line). Since the curves for the different temperatures fall onto 
a master curve, it is sufficient to focus on the curve for the lowest temperature in order to 
test whether the first scaling law holds. The results of our analysis are shown in Fig. 16 and 
17 for q = q max and q = q m in, respectively. Let us discuss q = q max first. As can be seen 
from Fig. 16, the von Schweidler law, with b = 0.55 = A = 0.76, the same value we used 
to fit the von Schweidler law to F s (q,t), fits the data over about 2.5 decades in time (long 
dashed line). Taking into account the next order correction (cf. Eq. 9) leads to a significant 
improvement of the fit for t > 4 • 10 4 (solid line). We have also fitted the numerical data 
with the /3-correlator with (dotted line) and without (short dashed line) correction. This 
was done by solving the equation for g_{t/t a ) (see Ref. 0) and making use of Eqs. (3) 
and (4). Here we encounter the same problem as we did for the fit with the von Schweidler 
law, because the /3-relaxation time t a can only be determined up to a T-independent factor. 
In order to proceed we have chosen for t a the position of the inflection point of F(q max ,t), 
indicated in Fig. 16 by a filled circle. The optimum values A = 0.76 and t a = 69 stemming 
from that fit were used for the similar fit of all the other correlators of the present paper 
and of Ref. [I5|, except for the orientational correlators with / = 1 ]T3. From the figure we 
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recognize that in the late /^-relaxation regime these fits are identical to the ones of the von 
Schweidler law, as it should be. In the early /3-relaxation regime, however, the /3-correlator 
fits the data much better than the von Schweidler law, since part of the approach to the 
plateau is fitted well also. Nevertheless, the critical law is not really observed, even if one 
includes the corrections to it [cf. Eq. (13)]. 

For q = q m i n it is not possible to obtain a good fit with the von Schweidler law alone, if 
h is kept fixed at 0.55. If b is used as a free parameter, a satisfactory fit is obtained but at 
the cost of a g-dependent von Schweidler exponent b q , in contradiction to MCT. Since the 
critical NEP obtained from this fit were rather structureless and did not qualitatively agree 
with the NEP obtained from solving the molecular MCT-equations P^] , we decided to keep 
b = 0.55 (obtained from 7^) fixed and to take the next order correction into account. With 
this approach we obtained good agreement between all NEP from our MD-simulation and 



those from the MCT-equations [P2|| . But even if the first correction to the von Schweidler 
law (cf. Eq. 9) is taken into account, the fit with Eq. 9 works well for only about one decade 
in time. On the other hand, the critical correlator (with A = 0.76) including the correction 
(on the von Schweidler side) fits the data over two decades. Nevertheless, this range is 
significantly smaller than the one found for q = q m ax which shows that corrections to the 
scaling laws might be more important for q min than for q max . 

The quantities f c (q), h(q) and h( 2 \q), depicted in Fig. 18 were determined in a similar 
way as for the self part, i.e. by keeping the value of b fixed to 0.55. Comparing these 
quantities with S(q) (Fig. 1) we find that f c (q) is in phase and both, h(q) and h^ 2 \q) in 
anti-phase with S(q). We note that the peak at q ~ 3 in f c (q) does not exist for S(q). 
However, whether this peak is real or just a statistical fluctuation can presently not be 
decided for sure. Since the static correlator S^q) has a pronounced maximum at q ~ 3 
(see Ref. [[TJJ), one might be tempted to relate this prepeak of f° to the translation- rotation 
coupling. 



The variation of f c (q) and h(q) with q resembles the one found, e.g., for hard spheres [E7 



for a binary liquid @,0] and that for water molecules . Interestingly the g-dependence 
of h^ 2 '(q), which is in phase with that of h(q), qualitatively agrees with that found by Fuchs 
et al. for hard spheres [^7]. However, in contrast to the hard sphere system, h^ 2 \q) does 
not change sign for the g-regime we have studied, but vanishes at q — q ma x at which it has a 
minimum. This fact explains why the (asymptotic) von Schweidler law fits the data rather 
well for q = q max , where h^ 2 '(q) ~ 0, but not for q = q m i n . To conclude we show in Fig. 19 
that also the variation of the a-relaxation time r q with q is in phase with S(q), which is in 



close analogy with, e.g., the hard sphere system [27 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have performed a detailed test of the predictions of ideal MCT (for 
simple liquids) for the translational degrees of freedom (TDOF) for a supercooled liquid of 
rigid diatomic molecules without head-tail symmetry. The Lennard- Jones interactions we 
used for the molecular system are not so different from those of the binary liquid |29] with 



80% A- and 20% B-atoms. The main difference between both systems of course is, that 
for the present case the 50% A- and 50% B-atoms are pairwise connected in order to form 
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diatomic molecules. Therefore the comparison of the dynamical behavior of both systems, 
allows to discuss the influence of the orientational degrees of freedom (ODOF) on TDOF. 

Our test was mainly concerned with (i) the existence of a single critical temperature T c 
and (ii) the validity of the two scaling laws of ideal MCT. This single critical temperature in- 
dicates a strong coupling between TDOF and ODOF. Although the mathematical structure 
of the MCT-equations for molecular liquids ||T9|-p^1 differs from the one of simple liquids, we 
expect that the two scaling laws hold for molecular liquids, too. This belief is based upon 
the fact that the first scaling law is always valid in a so-called type B transition (i.e. the 
NEP changes dis continuously at the transition), the type of transition that is relevant for 
structural glasses || . Furthermore the existence of a single transition temperature for linear 
and arbitrary molecules without additional symmetry can also be proved for the molecular 
MCT-equation g. 

Taking T c and 7 as free parameters, the a-relaxation times and r q for q = q max 
and q = q m in can be fitted with a power law. The resulting T c 's differ from T c = 0.475 
(deduced from D(T)) by less than 2%. The same value of T c was found for the orientational 
correlation functions in Ref. Jl6[ and we therefore conclude that our data is compatible with 
the existence of a single transition temperature. The sharp transition at T c , as predicted 
by ideal MCT, is however smeared out, due to ergodicity restoring processes which can 
be accounted for by the extended MCT [|]j]]. These processes are often associated with a 
hopping of the atoms or molecules, as it has been demonstrated for a binary liquid |[35|| . An 
evidence for such hopping processes is the occurrence of a second peak in G s (r,t) at r « 1 
for large times ||35|| . But no such peak is observed in our results, even beyond the r-range 
shown in Fig. 6. The absence of hopping processes with respect to the TDOF was also found 
in the MD-simulation of CKN ||. Therefore only orientational jumps remain, which indeed 
have been shown to be present |]16| . Comparing the different values of 7, one finds that 
those determined from F s (q,t) and F(q,t) fluctuate around 2.55. The same is more or less 
true for the purely orientational correlators Ci\t) and Ci(t) [17] for I > 1, but not for I = 1 

(s) 

where 71 = 1.66 and 71 = 1.52 was found. This result demonstrates that, contrary to the 
prediction of ideal MCT, 7 is not a "universal" exponent for our system. Such a conclusion 



was already drawn for binary liquids [29]. We note however, that this discrepancy does not 



seem to exist for the MD-simulation of water [ITS 



For all the correlators related to the TDOF, the second scaling law is reasonably well 
fulfilled. This is in contrast to our findings for c[ (t) and C\(t) ]TB[. Only for relatively 
large values of I do the Cj (t) show the second scaling law. 

Regarding the first scaling law we can say that it works well for F s (q, t) and F(q, t). For 
q = Qmin we have demonstrated that the next correction to the von Schweidler law becomes 
important. We stress that a consistent description of our data in the framework of MCT is 
only possible by taking this correction into account. If this is not done, a g-dependent von 
Schweidler exponent and a rather structureless g-variation of the orientational NEP (which 
will be discussed in the following paper ||18|| ) results. The critical law is not observed for any 
of the correlators, even if the next correction is taken into account. This is probably due to 
a strong interference of the critical dynamics with the microscopic motion. 

In summary, we conclude that, with respect to ideal MCT, the TDOF of our molecular 
system behave quite similar as found for a binary liquid [29]. Hence, the qualitative features 
of the dynamics of TDOF is not altered by its coupling to ODOF. Although a single transi- 
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does not fit into the framework of ideal MCT 
-reorient at ional jumps 



tion temperature for TDOF and ODOF can be specified, part of the orientational dynamics 
as measured, e.g. by C[ (t), 
this is only due to the 180 c 



161. Whether 



16], is presently not clear. 
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FIGURES 
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FIG. 2. S(q) (bold solid line) and the partial structure factors S a p versus q for T = 0.477; 
Saa (solid line), Sbb (dotted line) and Sab (dashed line). 
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FIG. 3. r-dependence of the atomic radial pair distribution function g a p(r) for T = 0.477; 
qaa (solid line), gBB (dotted line) and gAB (dashed line) 




FIG. 4. Radial distribution function g cm (r) for the center of mass positions and for T = 0.477 
(solid line), 0.63 (dotted line), 1.1 (short dashed line) and 2.0 (long dashed line). 
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FIG. 5. Temperature-dependence of the translational diffusion constant D; numerical data 
(squares, including error bars), power law (dotted line). The solid line serves as a guide for the 
eye. Inset: time dependence of the mean squared displacement for all temperatures investigated. 




FIG. 6. r-dependence of 4vrr 2 G s (r, t) for t = 0.4, 1.7, 5.88, 24.3, 101, 416, 1723, 5000, 2.5 • 10 4 , 
10 5 and 3.5 • 10 5 (from top to bottom) and for T = 0.477. 
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FIG. 7. (G s (r,t) - G s (r,t'))/(G s (r',t) - G s {r',t')) versus r for fixed r' = 0.105 and t' = 5000 
for 3.44 < t < 1723 and for T = 0.477. Inset: this ratio for the same r', t' and T, but for a 
much larger time span 0.4 < t < 3.5 • 10 5 . The times t are equidistant on a logarithmic time 
axisit ~ 0.40 • 2 n with n = 0, 1,2, . . .. Dotted lines: short times; solid lines: /3-relaxation regime; 
dashed lines: long times. 




FIG. 8. The non-Gaussian parameter 0:2 (t) versus time for all temperatures investigated. 
Inset: the maximum position t max versus T — T c with T c = 0.475. 
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FIG. 9. r-dependence of Gj, normalized by the average density p, for the same times as in 
Fig. 6 and for T = 0.477. 
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FIG. 11. (a) Time-dependence of the self part F a for q = q max = 6.5 for all temperatures 
investigated, (b) F s (q max ,t) versus rescaled time. The eight curves (solid lines) on the left refer to 
the eight lowest investigated temperatures and the one farthest to the right (dashed-dotted curve) 
corresponds to the highest temperature T = 5.0. The von Schweidler law including corrections (cf. 
Eq. 9) is represented by the dotted line, and a KWW-fit is shown dashed. 
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FIG. 12. a-relaxation times r q versus T — T c : open diamonds from F(q max , t), filled diamonds 
from F(q m i n ,t), open squares from F s (q max ,t), filled squares from F s (q m i n ,t). The dashed lines 
are a guide for the eye and both solid lines represent fits with a power law. 
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FIG. 15. F(q max ,t) (a) and F(q m i n ,t) (b) versus rescaled time. The eight curves (solid lines) 
on the left refer to the eight lowest temperatures and the one farthest to the right (dashed-dotted 
curve) corresponds to the highest temperature T = 5.0. The dashed line is a fit with a KWW-law. 
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FIG. 16. F(q max ,t) versus t for T = 0.477: numerical data (bold dots), /3-correlator with 
A=0.76 (short dashed line), /3-correlator with A= 0.76 including the next order corrections in the 
von Schweidler regime (dotted line), von Schwcidler law with b = 0.55 (=A = 0.755) (long dashed 
line), von Schweidler law with b = 0.55 including the next order corrections (solid line). 
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FIG. 17. F(q m i n ,t) versus t for T = 0.477: numerical data (bold dots), /3-correlator with 
A=0.76 (short dashed line), /^-correlator with A= 0.76 including the next order corrections in the 
von Schweidler regime (dotted line), von Schweidler law with b = 0.55 (=A = 0.755) including the 
next order corrections (solid line). 
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FIG. 18. Wave vector dependence of f c (q) (filled squares), h(q) (open circles) and h^ 2 \q) 
(open diamonds). 
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FIG. 19. Comparison of S(q) (solid line) for T = 0.477 and the a-relaxation time T q (filled 
squares) for T = 0.477, where T q has been multiplied by 10~ 5 ; inset: log S(q) and logr g versus q. 
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